Quantum— Classical Correspondence: Controlling Quantum Transport by State 

Synthesis in Ion Traps 

Juan F. Poyatos^ and Gonzalo Garcia de Polavieja^Q 

^ Centre for Quantum Computation, Clarendon Laboratory, Parks Road, Oxford 0X1 3PU, United Kingdom 
^Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford 0X1 3QC, United Kingdom 

(February 1, 2008) 



0^ 

On 

Oh- 
< 



> 
ON 

O 

o 

On 
On 

Oh! 

:3 



X 



A procedure to enhance the quantum-classical correspon- 
dence even in situations far from the classical limit is pro- 
posed. It is based on controlling the quantum transport be- 
tween classical regions using the capability to synthesize ar- 
bitrary motional states in ion traps. Quantum barriers and 
passages to transport can be created selecting the relevant 
frequencies. This technique is then applied to stabilize the 
quantum motion onto classical structures or alter the dynam- 
ical tunneling in nonintegrable systems. 
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Systems for which predictions can be made using clas- 
sical mechanics can show quantum mechanical properties 
under suitable experimental conditions. The observables 
measured can be of a non-classical class as in Bell experi- 
ments. States can also be chosen to deviate from classical 
behavior, an idea that we will exploit in the present Let- 
ter. Despite this, the correspondence principle continues 
to be a guide to the study of the interrelation between 
quantum and classical mechanics, specially applied to the 
study of atomic and molecular systems . More recently, 
a new class of experiments for which Hamiltonians can 
be engineered and detailed properties can be monitored 
have allowed to apply this correspondence in more detail. 
Cold atoms experiments have shown the possibility of 
inducing quantum dynamics with a particularly appeal- 
ing classical limit, the J-kicked rotor. These experiments 
observed dynamical localization, accelerator modes and 
only recently the effect of noise and dissipation ||]. In 
addition there have been theoretical proposals of exper- 
imental configurations in ion traps to study dynamical 
localization revivals Q and state sensitivity |^,^. 

In this Letter we show how to enhance the quantum- 
classical correspondence by controlling and monitoring 
quantum transport in an ion trap. The control is 
achieved not only by engineering the Hamiltonian but 
most importantly by state synthesis. To monitor the rel- 
evant effects we use tomography and simpler alternative 
techniques. Trapped ions advantages to both control and 
monitor have already been used for the study of other 
aspects of quantum systems, like the quantum Zeno ef- 
fect , possible nonlinear variants of quantum theory [p[ , 
reservoir engineering and quantum computation |l0| ]. 

To see the relative contribution of the classical back- 
bone and the purely quantum effects consider the Q dis- 



tribution that it is measured in tomography Its 
continuity equation is of the form 



Z{V,Q,t,n)=Q, 



(1) 



with Q{a;t) = being \a) a coherent state, 

V the potential and {,} the classical Poisson bracket. 
This continuity equation is the classical Liouville equa- 
tion plus a quantum term Z{V^Q^t^?i) given by Z = 
-|lm{(^|a)V (x/2 + ihd/dp) An ideal experi- 

mental setup to study the quantum-classical correspon- 
dence would then have the possibility to control all three 
dependencies of Z, i.e., the potential, the quantum state 
and an effective Planck constant. Moreover, it should be 
possible to have measurable quantities showing the rela- 
tive importance of classical and quantum contributions, 
e.g. the Q distribution itself. All of the above conditions 
can be fulfilled in the case of an harmonically trapped 
ion. This is possible due to the interaction between the 
internal (electronic) and external (vibrational) degrees of 
freedom of the ion by means of laser pulses in resonant 
and non-resonant regimes |l4| . 

In the quantum-classical boundary both the classical 
Liouville and the quantum terms in the continuity equa- 
tion in (|^) are relevant. The quantum contribution to the 
transport is state dependent and the final flow is diverted 
from the classical flow by an amount that depends on the 
quantum state. We first pick out the classical backbone 
structure and then study how to manipulate the quantum 
transport depending on the synthesized state. We need 
to construct a family of Hamiltonians that have regular, 
mixed or chaotic phase space. In the trapped ion setup 
we make use of the harmonic delta kicked Hamiltonian, 
first introduced in ||], describing a harmonic oscillator 
periodically perturbed by non-linear position dependent 
delta kicks. We consider a single trapped ion in a har- 
monic potential with two internal levels |e) and \g) 
with transition frequency coq interacting with a time de- 
pendent laser pulse of near-resonant light of frequency 
ul which is rapidly and periodically switched. For suf- 
ficiently large detuning S = uoq — uol the excited state 
amplitude can be adiabatically eliminated. The Hamil- 
tonian in this limit is given by [|| 

oo 

H = Ho + K[cos(2kLx) + l]\g){g\ ^ 6{t-nT), (2) 



1 



where Hq is the harmonic oscillator Hamiltonian, K 
is the kick strength related to the laser beam power, 
kL — 27r/A the laser wave number, and r the time be- 
tween kicks. A particularly important parameter is the 
so-called Lamb-Dicke parameter r] = kL^/fip^mv^ where 
V is the oscillator frequency and m the mass. Thus by 
varying r] it is possible to change the effective h of the 
system, i.e. doing it more classical. 

In the following the relevance of the classical backbone 
is shown by predicting the time averaged quantum Q 
function by means of classical theorems. We write the 
potential for a phase space region as V = + Vi, with 
Vo and Vi the unperturbed an perturbed potentials that 
in general have different phase space topologies. The 
time averaged dynamics of an initially synthesized mo- 
tional Fock state is shown to be predicted from knowl- 
edge of the classical solution for Vb and the Kolmogorov- 
Arnold-Moser (KAM) and the Poincare-Birkoff (PB) the- 
orems These two theorems in conjunction mean 
that, increasing the perturbation Fi, the phase space 
tori break in increasing irrationality values of the ratio 
of their winding numbers r = 001/002 with coi and C02 the 
two frequencies of a given torus jl^. When the winding 
number is sufficiently close to the rational number l/s 
the torus breaks into alternative ks stable and unstable 
points, with integers [|5|. 
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FIG. 1. (a) Stroboscopic classical map with k = 0.3. The 
rational number to which the winding number approaches is 1/4 
for the initial condition (x,p)=(0.22,0). (b) Qt function with 
k = 0.4 and ?7 = 0.25 (see text), (c) Q function averaged over 
three different tomographic measurements after applying 21, 22 
and 23 kicks. All parameters as before, (d) Stroboscopic classical 
map with k = 0.4. ur = 1.8 for all the cases, x and p are in 
units of A and miyX in all discussions and k = V2Kr]^ /h. 

Fig. 1(a) shows the classical stroboscopic map for Vq. 
Following the values of the winding number increasing 
the perturbation we determine which torus is breaking 



ffist. We thus predict the averaged quantum behavior 
with total potential V = Vq -\- Vi using only the classi- 
cal information from the values of the winding numbers 
with V = Vq and the classical theorems. Then a |n) = |7) 
Fock state is predicted to have an averaged Q function 
with four maxima (minima) on the classically stable (un- 
stable) points when increasing the perturbation of the 
system. Fig. 1(b) shows the time-averaged Q function 



Qt = lim 

M^oo Mr 



M-l 



J2QimT) = J2Ml^)?\{m\', (3) 



m=0 



with l/i) the Floquet states fuffilling t/(r)|/i) = e^^^^^l/i), 
where U{r) is the time evolution operator referring to 
one period r and ^ is a given initial state. The four 
predicted maxima are clear in this Figure. Fig. 1(c) dis- 
plays the Q function averaged at three consecutive times, 
as it could be measured experimentally using three to- 
mographic measurements showing the same struc- 
ture |1^. The classical stroboscopic map for the per- 
turbed case in Fig. 1(d) shows clearly that Qt reveals the 
classical backbone in many details. Any other Fock state 
prepared would scan in the same way different structures 
of this or alternative classical maps. 

We now propose to go a step further in the study 
of the quantum-classical correspondence controlling the 
relevance of the quantum and classical terms in ([|) by 
state synthesis. We want to prepare a state initially 
localized in a classical region A with a barrier or a 
passage for transport to a different classical region B. 
Take as a starting state a coherent state localized on 
^7 ^a(O). Its averaged transport to region 5, repre- 
sented by another state is given by P(|(/)/^), |(/>b)) = 
limM^oo 7v?7 5])m=o ^^^^ we can write 

as P{\ct^A).\M) = E;.K^BlM)l'K^A(0)|/i)P using the 
Floquet basis. For P(|(/)^), |(/>b)) not to be zero, the 
states (/>a(0) and (j)B niust have nonzero overlap with 
common Floquet states denoted as {|/i)AnB}- To form a 
new state that minimizes, (^^P, or maximizes, (p\^ ^ the 
transport to region B we eliminate (enhance) from ^a(O) 
the Floquet components {fiAns} as 

= N.ilMO)) ± E rMil^lMm, (4) 

with = or 1 when the corresponding weight p 
is smaller or greater than a value Ctoi, respectively and Ni 
will be normalization constants from now on. The value 
Ctoi is chosen to be the minimum possible subject to the 
condition that the new region is sufficiently close to 
the initial region of localization A. We propose then the 
synthesis of the states 



(5) 



n=0 
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with {\n)} the Fock basis and where we consider only an 
experimentahy feasible maximum number of Fock states 
nexp and amplitudes {n\ip^P) chosen to approximate 
the theoretical state (0). 
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FIG. 2. (a) Qt function for coherent state, (b) Qt 
function for modified state, (c) Upper Plot. Correlation 

function for 0^, (dot) and (dashed). Lower Plot. 

Autocorrelation function for state and 0a coherent state 
(dot), (d) Stroboscopic classical map. (e) Q for a initial coher- 
ent state centered at (xa^Pa) = (0? 0) after two kicks, (f) As 
before for a ^^'^ state. All figures obtained with the parameters: 
k = 2.7, z^T = 1.7 and 77 = .5. See text for more details of each 
case. 

The enhanced quantum-classical correspondence is 
thus obtained even far from the classical limit by choos- 
ing classical regions A and B and controlling the trans- 
port between them as described above. Some comments 
for this construction are necessary: (a) We can give up 
the condition that the initial state must be localized in 
region A' slightly bigger than the region A of the min- 
imum uncertainty wavepacket. For some experiments it 
is interesting to have a larger variance of the initial state 
or an additional localization in a different region A^\ (b) 
If we only want to modify the transport to region B and 
not to other regions C we have to use = in (||) for 
the Floquet states that have important overlap with a 
state 0c centered on C. (c) The states localized in the 
regions A, B and C do not need in general to be mini- 
mum uncertainty states but we have the freedom to use 
any states related to these phase space regions. 



We now discuss first a general example of modified 
transport properties of a given state and we then study 
stabilization, and modification of dynamical tunneling 
as two interesting examples. We want to synthesize a 
state ^^/^ initially localized on (xa^Pa) = (0.1,0.2) that 
avoids transport to {xb^Pb) = (0.3,0.3). We use co- 
herent states for (j)A and 0b and a value Ctoi = 0.001 
so all the Floquet states with important overlap with 
(j)B have = 1 in (^. An experimentally realizable 
state with nexp = 10 can be constructed with over- 
lap \{'^^P\(f^!^)\'^ = 0.9 to the theoretical state. De- 
spite both wavefunctions {a\(j)A) and being ini- 
tially localized around A, with = 0-^^, it 
is clear from Fig. 2(a) and (b) that their Qt functions 
are very different. The Qt function for ^^/^ shows a 
clear hole at the location B, in contrast to the min- 
imum uncertainty state 0/^. The correlation function 
P(mr) = |(0B|t/^(r)|c&)p is shown in Fig. 2(c) for 
$ = 0A and ^^P, the latter showing a significant de- 
crease. These correlation functions can be measured ex- 
perimentally by applying a displacement associated to 
(j)B and then measuring the Fock ground state population 
both steps possible to implement in a trap [|l3[ . Success- 
ful results can also be obtained for the state that 
presents a maximum on region B in Qt, see Fig. 2(c). 

Particular instances of the states are of special 

relevance. The stabilization of the vibration of the ion 
can be understood as a modified transport problem. In 
this case A is the region of localization and B the rest 
of the phase space and we are interested in construct- 
ing a state that we will name for the stabilization 
case as . An alternative way to understand expres- 
sion ^ for this situation now in terms of the dynam- 
ics is given in the following. The autocorrelation func- 
tion of a typical state |(0a(O)|0a(^))P will show an ini- 
tial maximum (however small) at > 0. We can now 
clean the state 0^ of the Floquet components that do not 
contribute to this maximum and therefore create a new 
state stabilized on A. Note first that a particular Flo- 
quet state l/io) can be obtained from the time-dependent 
vector |0a(^)) = |/i)(/i|0A(O)) exp(-ze^t) (a solution 
of the Schrodinger equation only for t = mr) as |/io) oc 

limT^oo G't,/.o with Gt,/.o = /Tt I^a(O) exp(ze^ot). 
A state related to the short term dynamics is then 
GtM.ujQ with uoq the value of uo that makes S{(jo) = 
/_!^^(it(0A(O)|0A(t))exp(zcjt) a maximum |l|]. This 
state GtM,uo can then be approximated as 

Ns^\^i){^i\4>Am, 



L^sta\ 

Wa' ) 



(6) 



with /ista the Floquet eigenfrequencies in the interval 
^^~Tt <^ < ^o + t^ with a weight |(/i|0A(O))p > agta- 
The value of agta is chosen maximum with the require- 
ment than the state (/p^f has a tolerable localization 
around A. The state in (||) synthesized to approxi- 
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mate (p^\f will then show an initial localization around A 
because it selects the short term dynamics and will recur 
continuously to A because it is made of very few selected 
Floquet states. In fact there are several maxima in S{(jo) 
and we can choose the value uoq with minimum number 
of Floquet states maximizing stabilization in this way. 

Using this stabilization procedure we have found en- 
hanced localization onto KAM tori, islands of size smaller 
than the effective fi^ cantori or unstable periodic orbits 
effects The following example of enhanced local- 

ization onto a classical unstable orbit can be realized 
experimentally. We consider the stroboscopic map in 
Fig. 2(d). Fig. 2(e) shows the Q function for a initial 
minimum uncertainty state, after two kicks centered 
at {xa^Pa) = (0,0). The state is spread from an un- 
stable periodic orbit along the unstable manifold. The 
stabilization is achieved in this case with a state cen- 
tered on the middle of the chaotic region at (0, 0) with 



0.72 and n, 



exp 



12, a single Floquet 



state. The stabilization achieved is clearly seen in the 
measurable autocorrelation function jl^], see Fig. 2(c) 
and Fig. 2(f) for the Q function of such stabilized state. 
This example constitutes a realizable experiment to di- 
rectly observe a quantum scar. |2Q] 




FIG. 3. (a) Stroboscopic classical map with k = 1.2 and 
ur = 27r/3. (b) Autocorrelation function for an initial coherent 
state, Floquet doublet (dot), and l^^,) state (dash). Parameters 
as before with 77 = 0.5. 

As a final example we show how to increase the dynam- 
ical tunneling associated to classically forbidden regions, 
see Fig. 3(a). An initial minimum uncertainty state lo- 
cated at (xa^Pa) = (—0.6,0) has dynamical tunneling 
and contributions from a unstable periodic orbit located 
at (xb-,Vb) = (—0.5,0.3). A state (H) eliminating the 
contributions of B is mainly composed of three Floquet 
doublets reflecting the oscillations due to quantum tun- 
neling. The autocorrelation function is shown in Fig. 
3(b) for a state (|) with {{(Pa^'^a^)? = 0-85. This can 
be measured by applying displacements associated to the 
initial regions of localization or more directly by inverting 
the unitary process which created the initial state. An 
example of stabilization in this case is choosing a state 
made of a single doublet where = 0.8. In 

this case the oscillations due to the tunneling are more 
clearly reflected. Both states with ngxp = 25 since their 
location in phase space implies higher Fock state contri- 



butions than previous cases. Choosing different doublets 
would give rise to different tunneling times. Finally no- 
tice that by increasing nexp, i-e. the overlap to the theo- 
retical state, all features discussed in previous examples 
will be more dramatically shown [|l^. 

In conclusion, we have presented a way to enhance the 
quantum-classical correspondence by modifying quan- 
tum transport in an ion trap. We have specially made 
use of the ability to synthesize arbitrary states of motion 
in a ion trap. 
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